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Abstract 

Comparisons of experimental data with theoretical predictions for collider processes containing 
hadronic jets rely on shower Monte Carlo event generators to include corrections to perturbative 
calculations from hadronization, parton showering, multiple parton collisions. We examine current 
treatments of these corrections and propose alternative methods to take into account nonperturba- 
tive effects and parton showering in the context of next-to-leading-order (NLO) event generators. 
We point out sizeable parton-showering corrections to jet transverse energy spectra at high rapidity, 
and discuss kinematic shifts in longitudinal momentum distributions from initial state showering 
in the case both of jet production and of heavy mass production at the Large Hadron Collider. 
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I. INTRODUCTION 



Phenomenological analyses of collider processes involving the production of hadronic jets 
rely on event simulation by parton shower Monte Carlo generators [l|, [2|. The subject of 
this paper concerns two different, common uses of shower Monte Carlo generators: one in 
which they are combined with hard scattering matrix elements via a matching scheme, e.g. 
at the next-to-leading order (NLO) 0, 0] in perturbative QCD; another in which they are 
used to obtain corrections to perturbative calculations due to hadronization, showering and 
multiple parton interactions (see e.g. js], |6[), with such correction factors then being applied 
to determine realistic predictions, which can be compared with experimental data. 

We begin in Sec. II by considering methods to evaluate the nonperturbative (NP) correc- 
tions to jet cross sections using shower event generators. We also estimate the corrections 
which arise from the initial-state and final-state parton showers and observe that they are 
sizeable (beyond NLO) in jet transverse energy spectra over the full range of rapidity. We 
propose a decomposition of the corrections to be applied to fixed NLO calculations, con- 
sisting of a truly NP contribution supplemented with a contribution coming from all order 
resummation via parton showers. 

Next, in Sec. Ill we investigate kinematic aspects of parton showers associated to combin- 
ing the approximation of collinear, on-shell partons with energy-momentum conservation. 
The main effect is an event-by-event shift in longitudinal momentum distributions whose 
size depends on the observable and on the phase space region, and increases with increas- 
ing rapidities. We illustrate this by numerical Monte Carlo results in different phase space 
regions for four specific examples of jet, heavy-quark, electroweak gauge-boson and Higgs 
boson production. First results on kinematic shifts have been presented in 0]. 

The approach of this work may be helpful to analyze corrections to finite-order perturba- 
tive calculations for jet observables from parton showering and nonperturbative dynamics. 
These encompass both final-state fragmentation effects and initial-state contributions asso- 
ciated with coUinearity approximations. Dynamical high-energy effects on jet final states, 
distinct from the ones discussed in this paper, have been emphasized in [8-10| due to non- 
coUinear contributions to parton branching processes. We note that both these results and 
the results in this paper stress the phenomenological relevance of more complete descriptions 
of QCD parton cascades in terms of transverse momentum dependent parton fragmentation 
and parton density functions lll-ll4|. Concluding comments on the results of this work are 
given in Sec. IV. 



II. MONTE CARLO NONPERTURBATIVE CORRECTION FACTORS 

In this section we consider methods to evaluate NP and parton shower correction factors. 
To be definite, we refer to the case of inclusive production of single jets at the LHC |l5[. In 
order to compare theory with experimental data corrected to stable particle level, Refs. s], |6[ 
supplement NLO perturbative calculations with NP corrections estimated from Monte Carlo 
event generators. Using leading-order Monte Carlo (LO-MC) generators [l|,|2|, the correction 
factors Kq are schematically obtained by [HI, [6| 

T^NP _ T^jips+mpi+had) ij.j{ps) x 
-"■0 ~ LO-MC I LO-MC ' \^) 

where {ps + mpi + had) and (ps) mean respectively a simulation including parton showers, 
multiparton interactions and hadronization, and a simulation including only parton showers 
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in addition to the LO hard process. Having only LO+PS event generators available, this 
is the most obvious way to estimate NP corrections to be applied to NLO parton level 
calculations. However, when these corrections are combined with NLO parton-level results, 
a potential inconsistency arises because the radiative correction from the first gluon emission 
is treated at different levels of accuracy in the two parts of the calculation. 

We here suggest that an alternative method which avoids this is to use NLO Monte Carlo 
(NLO-MC) generators to determine the correction. In this case one can consistently assign 
correction factors to be applied to NLO calculations. Moreover, this method allows one to 
study separately correction factors to the fixed-order calculation due to parton showering 
effects. To this end, we introduce the correction factors K^^ and K^^ as 

irNP _ j.j{ps+mpi+had) , i^jips) /r,\ 
^ NLO-MC I^^NLO~MC ' K'^ J 

L-PS _ t\t(ps) /yy(0) /o\ 

~ NLO~MCl NLO~MC ' \^ ) 

where the denominator in Eq. ([3]) is defined by switching off all components beyond NLO in 
the Monte Carlo simulation. The difference between the correction factors in Eqs. ([T]) and ([2]) 
comes primarily from the way in which the multiple parton interaction (MPI) contribution 
is matched to the NLO calculation. MPI processes have typical transverse momentum scales 
smaller than the scale of the hard process, which may be defined as the average transverse 
momentum of the hard partons. This however is different in LO and NLO calculations, 
giving rise to non-negligible numerical differences, which we will show below. The correction 
factor in Eq. ([3]), on the other hand, is new. It singles out contributions due to parton 
showering. This correction factor has not been considered in earlier analyses. We show 
below its numerical significance. We anticipate that taking properly into account these 
showering corrections can be relevant in fits for parton distribution functions using inclusive 
jet data. 

In Fig. [T] we compute results for the NP correction factors in Eqs. ([T]),([2I) to jet transverse 
momentum distributions. We define jets using the anti-k^ algorithm 16| with jet size 
R = 0.5 and R = 0.7. We plot the results versus the jet transverse momentum pt for different 
regions in the jet rapidity y. We show K^^ as obtained using the NLO event generator 
POWHEG [13] and compare it to the result obtained at leading order from Pythia ^ (tune 
Z2 [l8| and CTEQ6L1 pdfs [HI). The curves in Fig. [1] illustrate the differences coming from 
the definition of the hard process. 

In Figs. [2] and [3] we compute the corrections from parton shower K^^ as obtained from 
Eq. ([3]) as a function of the jet pt for different values of R and different rapidities y. Fig. [2] 
shows the contributions coming from initial state and final state parton shower separately. 
In general the effect from parton shower is largest at large \y\, where the initial state parton 
shower is mainly contributing at low pt, while the final state parton shower is contributing 
significantly over the whole pt range. 

While the NP corrections studied in Fig. [1] become vanishingly small at sufficiently large 
Pt, the showering correction in Figs. |2] and |3] gives finite effects also for large px ■ The 
dependence of this effect on the value of rapidity y will influence the shape of jet distributions 
and comparisons of theory predictions with experimental data. In particular, if the showering 
correction factor is not consistently taken into account, besides the NP corrections, this may 
affect the determination of parton distribution functions from data sets including jets. 

Note that in [s], [oj NP correction factors Ko are applied to the NLO calculation [20[ , and 
the data comparison shows that the NLO calculation agrees with data at central rapidities. 
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FIG. 1. The NP correction factors to jet transverse momentum distributions obtained from Eq. 
and Eq. using Pythia and POWHEG respectively, for \y\ < 0.5 and 2 < |y| < 2.5. Left: 
R = 0.5; Right: R = 0.7. 




FIG. 2. The initial and final state parton shower correction factor to jet transverse momentum 
distributions, obtained from Eq. (0j using PowHEG for \y\ < 0.5 and 2 < \y\ < 2.5. Left: R = 0.5; 
Right: R = 0.7. 
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FIG. 3. The parton shower correction factor to jet transverse momentum distributions, obtained 
from Eq. (0j using POWHEG for \y\ < 0.5 and 2 < |y| < 2.5. Left: R = 0.5; Right: R = 0.7. 



while increasing deviations are seen with increasin 

Pt ,A second comparison is performed in js 

tions 



rapidity at large transverse momentum 
with NLO-matched Powheg calcula- 
17|, showing large differences in the high rap^ity region between results obtained b;v 



tions [iY|j, snowing large diiierences m tne nign rapiditv region between results oDtamea by 
interfacing Powheg with different shower models [ll, Hj and different model tunes 18, 21|. Fl 
Motivated by this observation, in the next section we consider more closely the kinematics 
of the initial state parton shower at high rapidity. 



III. INITIAL STATE SHOWERING AND KINEMATIC SHIFTS 



Let us recall the physical picture [10[ of jet production at high rapidity (Fig. S]) based on 
QCD high-energy factorization 23|. Take the incoming momenta pi and p2 in Fig. HJin the 
plus and minus lightcone directions, and parameterize the exchanged momenta ki and ^2 in 
terms of purely transverse four-vectors k±i and k±2 and longitudinal momentum fractions 
Xi (collinear) and (anti-collinear) as ki = xipi + k^i + XiP2, k2 = X2P2 + ^±2 + ^2Pi- To 
single-logarithmic accuracy in the jet rapidity and the jet transverse momentum, we may 
approximate ki and ^2 using strong ordering in the longitudinal momenta, and get [lo[ 



ki ^ XiPi 



k. 



X2P2 



k 



±2 



Xi > X2 



(4) 



The physical picture corresponding to the factorization [10|, |23| consists of the scatterin g of 



a highly off-shell, low-x parton off a nearly on-shell, high-x parton. The calculations [10|, 122 



^ Further discussion of parton showering effects on high-rapidity jets may be found in 
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embody this picture through the longitudinal and transverse momentum dependences of both 
perturbative and nonperturbative components of the jet cross section, denoted respectively 
by a and $ in Fig. HI In what follows, however, we will not use the specific content of these 
calculations, but we will simply use the underlying physical picture as a guidance to examine 
kinematic effects of coUinear approximations. 




FIG. 4. Factorized structure of the jet cross section at high rapidity. 



In the light of this picture, let us consider the NLO-matched shower Monte Carlo calcu- 
lations, following In the Monte Carlo event generator first the hard subprocess events 
with full four-momentum assignments for the external lines are generated. In particular. 



the momenta kj^^ {j = 1,2) of the partons initiating the hard scatter are on shell, and are 
taken to be fully coUinear with the incoming state momenta pj, 

kf=x,p, (j = l,2) . (5) 

Next the showering algorithm is applied, and complete final states are generated including 
additional QCD radiation from the initial-state and final-state parton cascades. As a result 
of QCD showering, the momenta kj are no longer exactly coUinear, 

kj^XjPj (j = l,2) . (6) 

Their transverse momentum is to be compensated by a change in the kinematics of the hard 
scattering subprocess. By energy-momentum conservation, however, this implies a reshuf- 
fling, event by event, in the longitudinal momentum fractions Xj of the partons scattering 
off each other in the hard subprocess. The size of the shift in xj depends on the emitted 
transverse momenta. 

Let us now focus on jets measured in the rapidity range y < 2.5 ^ and examine the 
effect of the kinematical shift in the longitudinal momentum fractions. To this end we 
compute the distribution in Xj from Powheg before parton showering and after parton 
showering 0]. Fig. [S] shows the distribution for one of the Xj partons. We plot the result 
before showering (Powheg) and the results of successively including intrinsic kt, initial 
state parton shower and initial-|-final state parton shower. The results are obtained using 
the Pythia parton shower (tune Z2 1^ and CTEQ6L1 pdfs [3]). This does not include 



multiple parton interaction and hadronization effects. 
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FIG. 5. Distributions in the parton longitudinal momentum fraction x before (POWHEG) and 
after parton showering (POWHEG+PS), for inclusive jet production at different rapidities for jets 
with pt > 18 GeV obtained by the anti-kt jet algorithm \1 d / with R = 0.5. Shown is the effect of 
intrinsic kf, initial (IPS) and initial+final state (IFPS) parton shower. 



We see from Fig. |5]that the kinematical reshuffling in the longitudinal momentum fraction 
is negligible for central rapidities but becomes significant for y > 1.5. This effect character- 
izes the highly asymmetric parton kinematics, which becomes important for the first time 



at the LHC in significant regions of phase space [10|. Since the perturbative weight for each 
event is determined by the initial Powheg simulation, predictions of matched NLO-shower 
calculations for observables sensitive to this asymmetric region can be affected significantly 
by the kinematical shift as shown in Fig. El Similarly, since the momentum reshuffling is 
done after the evaluation of the parton distribution functions, the kinematical shift can af- 
fect predictions also through the pdfs. It will be of interest to examine the impact of this 
phase space region on total cross sections as well. 

Let us next consider the case of bottom-flavor jet production [24, 25|. The LHC measure- 
ments \24i . |25| are reasonably described by NLO-matched shower generators Mc@NLO 26 



and Powheg 23| at central rapidities, and they are below these predictions at large ra- 
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pidity and large px- In Fig. |6]we consider S-jets in different rapidity regions 2J] and plot 
the gluon x distribution from POWHEG before parton showering and after including various 
components of the parton shower generator, similarly to what is done above for Fig. [51 We 
use the Pythia parton shower (tune Z2 [18|], here including hadronization to identify the 
B-jet). We observe similar shift in longitudinal momentum with increasing rapidity as in 
the inclusive jet case. 





FIG. 6. Production of b-jets: distribution in the parton longitudinal momentum fraction x, before 
and after parton showering, for different rapidity regions. Shown is the effect of intrinsic kt , initial 
(IPS) and initial+ final state (IFPS) parton shower. 



In Fig. [7]we consider Drell-Yan (DY) production in the mass range 16 < m^jy < 166 GeV 
and perform a similar study to what is done above for jets. In this case too we find that the 
effects of the kinematical reshuffling in x evaluated from Powheg become non-negligible 
away from the central rapidity region. The double peak structure in Fig. [7] comes from the 
continuum DY production in addition to production. It will be of interest to investigate 
the kinematic reshuffling effect along with the forward Drell-Yan enhancements discussed 
in H 
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FIG. 7. Drell'Yan production with 16 < m/jy < 166 GeV: distribution in the parton longitudinal 
momentum fraction x before and after showering. Shown is the effect of intrinsic kt, initial (IPS) 
and initial+ final state (IFPS) parton shower. 



Finally we consider Higgs boson production in Fig. [H]for 110 < mmggs < 130 GeV. We 
observe a smaller effect at ^/s = 7 GeV than in the previous cases since the x-range is 
limited by the Higgs mass. 

Fig. |9] summarizes the results in Figs. [SHHl for the ratio of the cross section obtained by 
PoWHEG after inclusion of parton showering to the cross section before parton showering, 
plotted for different processes. In Fig. [10] we plot this ratio for Higgs boson production at 
different ^/s energies of 7, 14 and 33 GeV. 

The longitudinal momentum shifts from parton showering computed in this section mea- 
sure effects from QCD radiation beyond perturbative fixed-order calculations, and provide a 
significant contribution to the correction factors in Sec. II. They affect initial-state showers 
and need to be consistently taken into account in calculations which are used to determine 
parton density functions. The origin of the kinematical shifts lies with the approximation 
of collinearity 01 on the partonic states to which the branching algorithms describing show- 
ers are applied. Although for explicit calculations we have used a particular NLO-shower 
matching scheme (Powheg), the effect is common to any calculation matching NLO with 
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FIG. 8. Higgs boson production with 110 < mniggs < 130 GeV: distribution in the parton longitu- 
dinal momentum fraction x before and after showering. Shown is the effect of intrinsic kt, initial 
(IPS) and initial+final state (IFPS) parton shower. 



collinear showers. In calculations using integrated parton density functions the correction 
factors studied in this paper have to be applied after the evaluation of the cross section (and, 
as remarked on earlier, this may induce systematic inconsistencies if these corrections are 
not taken into account properly). On the other hand, this is avoided in approaches using 
transverse momentum dependent PDFs lll-ll4l. 1281] from the beginning (TMDs or uPDFs), 



as is done for example in the Cascade event generator [29 



IV. CONCLUSIONS 

Theoretical predictions for high-energy collider processes containing hadronic jets require 
supplementing finite-order perturbative calculations with parton showering and nonpertur- 
bative corrections. In this paper we have studied methods to treat parton showering and 
nonperturbative corrections in the context of matched NLO-shower event generators. 

We have pointed out potential inconsistencies in current approaches which on the one 
hand apply NP correction factors from leading-order Monte Carlo generators to NLO parton- 
level predictions and on the other hand fail to include showering corrections. We have 
proposed methods to address these deficiencies by using consistently available NLO Monte 
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FIG. 9. Ratio of the cross sections obtained with POWHEG after and before inclusion of initial + 
final state parton shower and intrinsic kt for the different processes. 



Carlo tools. We have shown that the differences in the predictions for jet cross sections 
induced by the modified approach we propose are significant in regions of phase space which 
are explored with hard probes for the first time at the LHC In particular, the nonpertur- 
bative correction factor K^^ introduced in Sec. II gives non-negligible differences at low 
to intermediate jet pt, and the showering correction factor K^^ of Sec. II gives significant 
effects over the whole px range and is largest at large jet rapidities y. 

Because of this y and pt dependence, taking properly into account NP and showering cor- 
rection factors changes the shape of jet distributions and affects significantly the comparison 
of theory predictions with experimental data. The numerical results we have presented show 
effects as large as 50 percent in regions of y and pt phase space relevant to jet measurements 
at the LHC. The showering correction factor K^^ in particular can affect the determination 
of parton distribution functions from fits to experimental data sets comprising inclusive jet 
measurements. 

We have investigated in closer detail the sources of the showering correction from initial- 
state and final-state effects. We have observed that the main initial-state showering effect 
comes from kinematical shifts in longitudinal momentum distributions 0] due to combin- 
ing coUinearity approximations with the Monte Carlo implementation of energy-momentum 
conservation constraints. We have examined the longitudinal shifts for specific processes in 
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FIG. 10. Ratio of the cross sections obtained with Powheg after and before inclusion of initial 
+ final state parton shower and intrinsic kt for Higgs production at different energies: y/s = 
7,14,33 TeV. 



Sec. III. This effect is largest for inclusive jets and 6-fiavor jets at the LHC in the higher 
rapidity bins. We have extended the study of longitudinal shifts to the case of Drell-Yan 
pair production by analyzing the Drell-Yan mass region 16 < mnY < 166 GeV and found 
that the shifts are non-negligible for Drell-Yan production at forward rapidities y > 2. We 
have also examined the case of Higgs boson production for 110 < mmggs < 130 GeV and 
found that the shifts are non-negligible at large rapidities at y/s = 7 GeV, and become more 
and more important at higher centre-of-mass energies. 

It will be interesting to study the impact of the effects discussed in this work on phe- 
nomenological analyses of LHC final states involving hadronic jets. We expect these effects 
to also influence determinations of parton distributions. Longitudinal momentum shifts 
can be avoided in formulations that keep track of non-collinear (i.e., transverse and/or 
anti-coUinear) momentum components from the beginning using unintegrated initial-state 



distributions |12l . |13| , also at parton shower level |29l . |30| . It will be interesting to investi- 
gate to what extent this can be exploited to construct approaches in which nonperturbative 
contributions such as multiple parton interactions, finite transverse momenta, hadronization 
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are consistently incorporated in parton branching event generators. 
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